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Abstract 

Melt-mixing in twin-screw extruders is a key process in the development of polymer composites. 
Quantifying the mixing performance of kneading elements based on their internal physical processes 
is a challenging problem. We discuss melt-mixing by novel kneading elements called "pitched-tip 
kneading disk (ptKD)". The disk-stagger angle and tip angle are the main geometric parameters 
of the ptKDs. We investigated four typical arrangements of the ptKDs, which are forward and 
backward disk-staggers combined with forward and backward tips. Numerical simulations under a 
certain feed rate and screw revolution speed were performed, and the mixing process was investigated 
using Lagrangian statistics. It was found that the four types had different mixing characteristics, and 
their mixing processes were explained by the coupling effect of drag flow with the disk staggering 
and pitched-tip and pressure flows, which are controlled by operational conditions. The use of a 
pitched-tip effectively to controls the balance of the pressurization and mixing ability. 

Keywords: Polymer Processing, Twin-screw extruder, Numerical simulation, Dispersive mixing, Dis- 
tributive mixing 

1 Introduction 

A twin-screw extruder (TSE) is widely used in various industries, including polymer processing, food 
processing, rubber compounding, pharmaceutical development, and processing other highly viscous ma- 
terials because multiple processes, for instance, heating, melting, compounding, blending, reaction, and 
devolatilization, can be conducted continuously in one operation. In the plastic industry, a TSE is used 
to develop products with certain desired properties by compounding and blending different types of poly- 
mers, fillers, and other additives. In this process, controlling the mixing in the TSE is critical to optimize 
and improve the productivity and functionality of products. 

The control of the mixing process is based on the design of the kneading screw element. In the past, 
among different kneading screws, kneading disks (KDs) were the most commonly used for general purpose 
compounding due to its mixing efficiency. One main geometric parameter of a KD is the disk-stagger 
angle that controls both the pumping capability and inter-disk leakage flow simultaneously. However, 
these two factors have different effects on the mixing process. 

To control these two factors separately, we have proposed a nov el design of a kneading screw element, 
called "pitched-tip kneading disk (ptKD)" ( Shigeishi et al. . 2006lh In the ptKDs, disk tips are pitched 



to the screw axes in contrast to the tips that are pitched parallel to the screw axes in conventional KDs. 
By arranging the pitched-tip angle along with the disk-stagger angle, the pumping capability and the 
inter-disk leakage flow are controlled independently. This flexibility makes it possible to better control 
distributive mixing and dispersive mixing than with conventional KDs. 
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For a specific application of novel screw or barrel designs, it is difficult to optimize extrusion conditions 
because there are a large number of parameters including geometric and operational parameters of the 
device and rheological parameters of the materials, which affect the flow pattern and mixing process. It 
is difficult to experimentally obtain information on the flow pattern and mixing process in a TSE because 
TSEs are operated under severe conditions, including high temperature and high pressure; moreover, the 
channel in TSEs in which materials are transported is highly complicated. Numerical simulation can aid 
in assessing material behaviors in a TSE. 

Thre e-dimensional fluid dynamic simulations of melt-mixing in TSE hav e been performed by several 



authors (IKaiiwara et al.1.119 96: Yoshinag a et all 12 000: Ishikawa et alll20 00a.b. 2001, 2002; Funatsu et al 



2002; Yang and Manas-Zloczowerl 1992c ICheng and Manas-Zloczowerl 11998c lYao and Manas-Zloczower , 
1998c Cheng and Manas-Zloczowerl ll997[|Bravo et all 120001 12004c Ijaffer et all l2000t lLawal and Kalvon , 



1995al lbl: iKalvon and Malikl , 12007c iMalik and Kalvoni [2005t lAlsteens et all 12004) . It is essential to per 



form three-dimensional simulations to observe the material flow and the mixing process in TSE because 
the process is mainly dominated by the three-dimensional structure of the channel. 

In this article, we investigate the flow characteristics and mixing process in melt-mixing with several 
types of ptKDs. We conducted three-dimensional numerical simulations of melt-mixing zone in a TSE. 
Material kinetics was solved with the Lagrangian method, and the characteristics of distributive and 
dispersive mixing in melt-mixing with several types of ptKD was investigated. 



2 Pitched-tip kneading disk 

One element of a ptKD is composed of several kneading blocks, like the conventional KD, but each 
kneading block of a ptKD has non-parallel tips to the screw axes, in contrast to conventional KDs. When a 
tip is pitched to pump fluid forward (backward) along the screw rotation, it is called a "forward (backward) 
tip" . When the disk-stagger angle is arranged to pump fluid forward (backward) along the screw rotation, 
it is called a "forward (backward) stagger" . According to these two main geometric characteristics, ptKDs 
are typically classified as one of four types: backward stagger and forward tip (Bs-Ft), forward stagger 
and forward tip (Fs-Ft), backward stagger and backward tip (Bs-Bt), and forward stagger and backward 
tip (Fs-Bt). Figure Q] shows the top views of the four typical types of ptKDs. 

The four types of ptKDs in Fig. Q] are composed of five blocks of one element and are chosen for our 
investigation. Based on the barrel diameter D — 69.0 mm, each self-wiping block (jBoov . 19781 ) had a 



flight radius of 0.481? and a width of 0.3D. Disk-stagger angles were set to ±40.685° for Fs-Ft and Bs-Bt, 
±49.315° for Fs-Bt and Bs-Ft, respectively. The leads of four types of ptKDs were set to ±6D. 
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Figure 1: Top view of typical pitched-tip kneading disks: (a) backward stagger and forward tip (Bs-Ft), 
(b) forward stagger and forward tip (Fs-Ft), (c) backward stagger and backward tip (Bs-Bt), (d) forward 
stagger and backward tip (Fs-Bt). 

3 Numerical simulation 

3.1 Basic equations 

Our simulation focuses on the situation in which the material is fully filled in the melt-mixing zone in 
a TSE. The flow of the polymer melt in the TSE was assumed to be incompressible, and the Reynolds 
number was assumed to be much less than unity so that the inertial effects could be neglected. We con- 
sidered the pseudo-steady state of mass and heat transfers under which the momentum and temperature 
fields instantly reach the steady state for a given channel geometry, and the time evolution of the fluid is 
totally governed by the screw rotation. With these assumptions, the governing equations become 



V • v = 0, 

= -Vp + V • r, 
pC p v ■ VT = kV 2 T + t : D, 



(1) 
(2) 
(3) 
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where p and t are the hydrostatic pressure and the deviatoric stress, respectively; p, C p , T and k are the 
mass density, the specific heat capacity, the temperature and the thermal conductivity, respectively; and 
D = Vt> + (Vu) T /2 is the strain-rate of the velocity field v, where (.) T indicates the transpose. 

3.2 Working fluid and constitutive equations 

We employ a homogeneous fluid as the working fluid because we are focused on characterizing the melt- 
mixing of different screw geometries. The working fluid is assumed to be a viscous shear-thinning fluid 
that follows the Cross-exponential model: 

r = 2 V D, (4) 
,( 7 ,T)= >»FoWT,T ) 

IKh ' i + ^To^^Toh) 1 -"' V ' 

H{T,T ) = exp[-l3(T-T )}, (6) 



7 = V2D : D, (7) 

The parameters of the Cross-exponentia l model were determin ed by fitting function to the experimental 
shear viscosity of a polypropylene melt ( Ishikawa et all [2000al) . as shown in Fig (2] 



The mass density, specific heat capacity, and thermal conductivity were assumed to be temperature- 
independent in the simulation: p = 735.0 kg/m 3 , C p = 2100.0 J/(kg-K), and k — 0.15 W/mK, respec- 
tively. 

3.2.1 Boundary conditions 

The no-slip condition on the velocity at the barrel and screw surfaces is assumed. The velocities at 
inlet and outlet boundaries were obtained according to the same cross-sectional velocity profile under the 
given volumetric flow rate of Q = 60 cm 3 /s («159 kg/h). The temperatures on the barrel surface and at 
the inlet boundary were set to 473 K and 453 K, respectively. The natural boundary conditions for the 
temperature equation in the exit boundary plane and the screw surface were assigned. 

3.2.2 Penalty method for the incompressibility condition 

To incorporate the incompressibility condition ([1} in momentum conservation ^ , the penalty method is 
utilized: the pressure in Eq.([5]) is eliminated by imposing the relation 

p = -Ary V • v, (8) 



with a constant A ^> l (|ZienkiewiczUl977h . The penalty term in the momentum conservation effectively 



acts as a very large bulk viscous stress; therefore, the resulting velocity field satisfies the incompressibility 
condition (TT]) in the limit as A — > oo. 

3.2.3 Spatial discretization and time evolution 

The channel in a TSE is decomposed by the hexahedral element. The velocity and the temperature are 
approximated by tri-quadratic interpolation with the 27-node in each hexahedral element. The momen- 
tum equations @ and © were discretized by the Galerkin finite element method. The temperature 
equation flSt was disc r etized by streamline- upwind/Petrov-Galerkin method ( Brooks and Hughesl 1982 : 



Marchal and Crochet . 19871 ) to stabilize the advection term. 



The time evolution of the velocity and temperature fields was constructed with the converged fields 
for every three degrees of screw rotation. This computation was the most costly. For reference, one 
simulation takes about 68 hours on a computer with an Intel® Core™ 2 Extreme processor. 
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Figure 2: Shear viscosity of a polypropylene melt dshikawa et all l2000ah . The solid lines are the Cross- 
exponential model ©-© fit. The parameters are T = 473 K, A(T ) = 1.70249 s, 77o(T ) = 32783 Pa-s, 
ra(T ) = 0.330793 (-), /3 = 0.0208143 K _1 . 
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3.2.4 Lagrangian point tracer 



The trajectories of the passive tracers were determined based on the solved velocity field. The time 
evolution of the tracer density was utilized to assess the flow pattern in the melt-mixing zone with a 
ptKD . The set of tracer trajectories was utilized to compute the residence time distribution ( LevenspielL 



19981 ) and to investigate the flow history and mixing process. The Lagrangian-history average of a 



quantity / over the trajectory of crth tracer is defined as 



—T a 1 ' T « 



fa " =JT j a dxS (x-X a ( S ))f( X , S ), (9) 

where T a and X a (.) are the residence time and position of the ath tracer, respectively, and 5(.) is the 

— T a 

Dirac delta distribution. The statistical distribution of f a characterizes the process in the melt-mixing 
zone. 

Initially, 2000 points were uniformly distributed in a certain section relative to the axial position, 
which is arbitrarily set in the second disk. They were advected until they reach the outlet section. When 
computing the tracer advection, some tracers went out of bounds because the time resolution of the 
velocity field was limited. To circumvent the effect of the lost tracers on the statistics, we set the number 
of initial points of 2000, which ensured that a sufficient number of points reached the outlet. In our 
advection computation, at least 93% of the tracers reached the outlet. 



4 Results and discussion 

4.1 Basic characteristics of ptKDs 

4.1.1 Drag ability of ptKD 

To assess the effect of pitched tips on the pressurization capability, the pressure drop from the inlet to 
the outlet was calculated under a fixed screw rotation speeds and various volumetric flow rates (Fig|3]). 
For a fixed volumetric flow rate, the pressure drop was largest in Bs-Bt, followed by Bs-Ft, Fs-Bt, and 
Fs-Ft. The results in Figj3]show that the disk stagger mainly controls the drag ability of ptKD and that 
the pitched tips only experience slight modification. It is noted that the pressure drop in Fs-Ft ptKD 
was negative for a certain range of the volumetric flow rate, which indicates that the pressure-driven flow 
on average occurred in the direction from the outlet to the inlet. 

4.1.2 Temperature profile 

Figure [4] shows the average temperature profile along the extrusion direction. For all types of ptKDs, the 
temperature rose from the inlet temperature to the barrel temperature over one block. Next, the average 
melt temperature increased two or three degrees beyond the barrel temperature due to viscous heating 
effects. Because the difference in the results with different types of ptKD was very small, we assume that 
the temperature did not have a significant effect on the mixing ability. 

4.1.3 Residence time and Lagrangian history of exerted stress 

FiguresOandOshow the joint probability density function (PDF) of the residence time and the Lagrangian- 
history average of stress magnitude under Q = 60 cm 3 /s and 200 rpm. To measure the stress magnitude, 
a second-order invariant of r is considered 0- Although r = 2i]D is traceless by the incompressibil- 
ity condition ([1]), a small but finite trr is unavoidable in numerical simulations with finite arithmetic. 

1 Note that although one of the authors (T.K.) denned "characteristics shear stress" as 

\J 4 (jxy + T yz + T zx) + ( r xx - T yy ) 2 + (r yy - T zz f + (j zz - T xx f , 

lwara ct al., 1996) which was used in Yoshinaga ct al ] 1120001) and llshikawa et al ] d200ll . l2002h . it is not an invariant of 
a second-rank tensor; therefore, it is not adequate to quantify the magnitude of the stress tensor. 
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Figure 3: Pressure drop averaged over screw rotation with a screw rotation speed of 200 rpm as a function 
of the volumetric flow rate. 
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Figure 4: Average temperature profile along the axial direction with a volumetric flow rate of 60 cm 3 /s 
and a screw rotation speed of 200 rpm. 
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Therefore, we use the traceless part of r to compute the stress magnitude, a, as 



-T> : T>, (10) 
T / = T-i(trr)J, (11) 



where J is the unit tensor. 



Counter-pump combinations: Fs-Bt and Bs-Ft In Fs-Bt and Bs-Ft, the pitched-tip angle and 
disk-stagger angle are in opposite directions, which we call the counter-pump combination. In the counter- 
pump combination, the joint PDFs of T a and ~da Ta are almost unimodal (FigJS|). These ptKDs are 
characterized by a broad residence time distribution and a rather small fluctuation of Oa Tc " . A small 
fluctuation of o r ^ Ta indicates that the average stress history of each trajectory is almost homogeneous; 
thus, homogeneous dispersive mixing can be expected for Fs-Bt and Bs-Ft. 



Co-pump combinations: Fs-Ft and Bs-Bt In Fs-Ft and Bs-Bt, the pitched-tip angle and disk- 
stagger angle are in the same direction, which we call the co-pump combination. In the co-pump com- 
bination, the joint PDFs of T a and Oa Ta are spread over the plane (Fig|5]), which indicates that the 
residence time and stress history of each tracer are highly inhomogeneous. 

For the Fs-Ft (FigJS](a)), the main peak indicates that large population of the tracer has a shorter- 
residence time and smaller average stress, which indicate that pipeline flow occurs where tracers extrude 
under fast velocity and a small strain-rate and mixing with surrounding fluid rarely occur. The other 
trajectories have larger average stress during longer residence time. 

For the Bs-Bt (FigJHKb)), the joint PDF of T a and Oa Ta is bimodal; although, the residence time 
fluctuation is smaller than those of other ptKDs. The Lagrangian-history average of the stress magnitude 
has much greater fluctuation than other ptKDs. The bimodal structure indicates that pipeline flow also 
occurs in Bs-Bt, but the group velocities of the pipeline and its surroundings are similar. Note that the 
indication of the pipeline flow cannot be observed based only on the residence time distribution. The 
joint PDF reveals the highly inhomogeneous tracer history in the co-pump combination. 

4.2 Transport mechanism and characterization of mixing capability 
4.2.1 Time evolution of tracer density 

Spatial density of the tracer as a function of time, p(x,t), is computed from the trajectories of the 
tracers, Xi(t) (i = 1, . . . ,N). To investigate the axial tracer kinetics, the tracer density is projected in 
the axial (z-) direction as 



1 N 

p(x,t)= lim ^(x-Xiit)), (12) 



i=l 



P( z ,t) = J dxdyp(x,t). (13) 

Figure l4.2.1l shows the evolution of the axial distribution for the four ptKDs, provides detailed information 
on the transport mechanism, and explains the stress- residence-time j oint PDF in the previous section. 



RTD has been discussed often in connection with the axial mixing (Ishi kawa et al. . 2002; Z hang et al 



20091) . However, a detailed investigation of the axial transport for ptKDs revealed that the source of the 



residence time fluctuation differs among the four types of ptKDs, and that broadness of RTD does not 
necessarily means good mixing capability. 
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Figure 5: Joint probability density of the residence time and Lagrangian-history average of the stress 
magnitude of Ea. (fT0|) with a volumetric flow rate of 60 cm 3 /s and a screw rotation speed of 200 rpm: (a) 
Fs-Bt ptKD, and (b) Bs-Ft ptKD. 10 
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Figure 6: The same plot as FigEJbut for (a) Fs-Ft ptKD and (b) Bs-Bt ptKD. 
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Figure 7: Time evolution of the tracer density on the axial direction for (a) Fs-Bt, (b) Fs-Ft, (c) Bs-Ft, 
and (d) Bs-Bt. The initial axial position of the tracer is depicted by the vertical line on the second disk. 
The different times correspond to the different lines indicated in the legend. 



Counter-pump combinations: Fs-Bt and Bs-Ft For the counter-pump combination, the inter-disk 
leakage flow in front of the pitched-tips enhanced the axial mixing. As for Fs-Bt in FigfTJa), a portion of 
the tracer was pumped backward at the inter-disk gap. A close investigation of each tracer's trajectory 
revealed that the compression flow in front of the backward tips caused this inter-disk backward pump 
and enhanced the axial mixing. The effect of backward pitched-tips and circumferential mixing at the 
intermeshing region promoted gross mixing over a whole channel of melt-mixing zone. 

As for Bs-Ft in FigjT^c), a portion of the tracer was pumped forward at the inter-disk gap by the 
forward pitched-tips. This effect also caused axial mixing as in the case of Fs-Bt. The broad RTDs 
observed in FigJS] are the result of axial mixing, and the unimodal stress-history is a result of gross 
mixing. 



Co-pump combinations: Fs-Ft and Bs-Bt For Fs-Ft in FigJTJb), it is observed that the unimodal 
axial density immediately split into a bimodal density. These two groups evolved with different axial 
group velocities. A close investigation of each tracer's trajectory revealed that the faster axial transport 
occurred at the roots of the kneading disks where it existed larger space to enhance the pressure-driven 
flow. The faster group had a pipeline flow characteristic in which mixing with surrounding fluid rarely 
occurred. Because the forward-pump ability of Fs-Ft ptKD is large, a negative pressure drop is generated 
in the volumetric flow rate (FigJ3]). The gross axial transport is a result of the competing drag and 
pressure gradient. The broader RTD in Fs-Ft (FigJSJa)) is a result of pipeline flow and the surroundings. 
Because the history of exerted stress was different between the pipeline region and the surroundings, 
7j^ T ° fluctuated significantly, as observed in Fig[6ja). 

As for Bs-Bt in FigJTJd), splitting of the axial density was also observed as in Fs-Ft, but the group 
axial velocities were closer. A close investigation of each tracer's trajectory revealed that one group had 
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a character of pipeline flow as in Fs-Ft. The large fluctuation of ~da Ta ( Fig|5{b)) was a result of pipeline 
flow and the surroundings. However, similar axial group velocities for the pipeline and its surroundings 
resulted in a narrower RTD (Figl^b)). 

4.2.2 Mixing capability: strain-rate type and maximum line-stretching rate 

For effective mixing, two points that are close initially must separate rapidly over time to promote 
points distribution in space, stretch a line clement, and stretch an area element. To quantify the mixing 
capability, we considered two aspects of this kinetic process: the potential line-stretching capability and 
the type of strain-rate field in which the line-stretching occurs. 
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Figure 8: Probability density of the Lagrangian-history average of the maximum line-stretching rate with 
a the volumetric flow rate of 60 cm 3 /s and a screw rotation speed of 200 rpm. Arrows on the abscissa 
indicate the mean values. 



To quantify the potential line-stretching capability, we used the maximum line-stretching rate 

\M(x,t)— max D : nn. 
n,\n\=i 



(14) 



The details of the quantity are described in [5] The statistical distribution of the Lagrangian-history 

T a 

average of the maximum line-stretching rate, (Am) q > is depicted in Fig|S]For co-pump ptKDs, namely 

T a T a 

Fs-Ft and Bs-Bt, PDFs of (Am)„ are broad. In particular, for Fs-Ft, a shoulder at large (Am) q 

T a 

is observed. The large fluctuation of (Am) q is a result of the pipeline flow. In the trajectory by 
the pipeline flow, (Ajvf) is small, and the larger (Am) q comes from the surroundings. The large 



fluctuation of (A 



Mi 



indicates that gross mixing might be inhomogeneous in co-pump ptKDs. It is 
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noted that the fingerprint of the pipeline flow is clearer in the PDF of the maximum line-stretching rate 
than in that of the residence time. 

T Q 

For counter-pump ptKDs, namely Fs-Bt and Bs-Ft, PDFs of (Am) q are unimodal, which indicates 
that no special structure exists in the flow pattern in counter-pump ptKDs, and mixing over whole channel 
is expected. 

T a 

It is noted that the mean values of (Am)„ indicated in Fig J5] are determined by the disk-stagger 
direction and are useless for characterizing the mixing capability. To characterize the mixing capability, 

T Q 

the fluctuation of (Am) q should be taken into account. 
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Figure 9: Probability density of the Lagrangian-history average of the strain-rate-type identifier with a 
volumetric flow rate of 60 cm 3 /s and a screw rotation speed of 200 rpm. 

To identify a type of strain-rate field, the strain-rate-type identifier as a function of the invariants of 
the strain rate, D(x,t), is introduced as 

max«, D : nn 

MX > t] = VD.D ' (15) 
where n is a unit vector. The strain-rate-type identifier takes values specific to the strain-rate state in 
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incompressible flow, 



(uniaxial elongational flow) 
Jd = \ \/l/2 (planar shear flow) (16) 
(biaxial elongational flow) 




The details of the strain-rate- type identifier are summarized in O 

The statistical distribution of the Lagrangian- history average of the strain-rate- type identifier, (/d)_ , 

T a 

is depicted in FigJT5] The modes of the PDFs for the four ptKDs are located near (/d) q = y 1/2 sa 
0.707, which is a general characteristic of screw extruders in which the contribution from circumferential 
shear flow is almost fjj = y^/2. A close look at the modes in the inset of Fig. |H] reveals that the 
ptKD characteristics were found. For co-pump ptKDs, namely Fs-Ft and Bs-Bt, the modes of their 
PDFs indicate that the history of the strain-rate type is planar shear flow. In contrast, for counter-pump 
ptKDs, namely Fs-Bt and Bs-Ft, the contribution of biaxial elongational flow at the pitched tips causes 

T a 

mode to shift to a lower value of (/£>)„ • The result suggests that the counter-pump ptKDs to generate 
the three-dimensional deformation, which promotes efficient mixing, are more effectively than co-pump 
ptKDs. 



5 Concluding remarks 

For melt-mixing with a novel mixing element, the "pitched-tip kneading disk (ptKD)", its basic charac- 
teristics and mixing capability were investigated based on numerical simulation of three-dimensional flow 
and the statistics of the Lagrangian tracer. The ptKD has two main geometric parameters: the pitched- 
tip angle and the disk-stagger angle. Four typical combinations of the pitched-tip angle and disk-stagger 
angle were used in the current research. 

The pressurization capability is mainly determined by the disk-stagger angle, and the pitched-tip 
provides minor modification. 

Flow patterns and transport mechanisms were studied according to the tracer density evolution. In 
co-pump ptKDs (Fs-Ft and Bs-Bt), pipeline flow exists that causes a large fluctuation in the Lagrangian 
history of several quantities, including the exerted stress and the line-stretching rate. In co-pump ptKDs, 
the pitched-tip enhances the pressurization and causes the axial transport to be inhomogeneous by gen- 
erating pipeline flow. 

In counter-pump ptKDs (Fs-Bt and Bs-Ft), the biaxial elongational flow occurs in front of the pitched- 
tip and enhances mixing over the channel. The result suggests that the main advantage of Fs-Bt and Bs-Ft 
is gross distributive mixing. It is expected that, the optimal balance of distributive mixing and dispersive 
mixing are found for various extrusion processes by adjusting the combinations of the pitched-tip angle 
and the disk-stagger angle in Fs-Bt and Bs-Ft. 

To quantify of mixing capability, the maximum line-stretching rate and the strain-rate-type identifier 
were introduced. Information on the flow patterns and mixing capability is reflected in the statistics 
of Lagrangian-history averages of the maximum line-stretching rate and the strain-rate-type identifier. 
These quantities are useful for analyzing the results of three-dimensional flow simulations. 



Other characterizations of mixing in three-dimensional flow u sing the Lyapunov exponent (|Lawal and Kalvon . 



1995a ) and intensity of segregation ( Lawal and Kalvon . 1995bl) have been suggested. It would be inter 



esting to compute these measures in ptKDs and to compare the results with those of other kneading 
elements. We are currently working on this project. 
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A Local line-stretching rate and type of strain-rate field 



The mixing ability of a flow field can be quantified by the line-stretching rate of an infinitesimally short 
line element I. The time evolution of / obeys for a short time period 



dl 

dl 



= 1 D 1 



(17) 



The stretching rate of the norm \l\ is derived from Eg. (TIT)) as 



i <yiz 

III dt 



I dl 

\f\ ' dt 

— n — 

\i\ ' ' izi 



(18) 



When the line-stretching rate D : 11/ \l\ is positive, an infinitesimally small distance nearby two points 
separate in exponential manner. 

The stretching rate of |Z| (p~8|) depends on the direction of I. To characterize the mixing ability of a 
given flow field, the maximum of Eq. (|18p or equivalently the maximum line-stretching rate, Am in Eq. (|14[) 
can be used. By definition, Am coincides with the largest eigenvalue of D. 

In the theory of dynamical systems in a bounded system, the average stretching rate based on the 
director advected by Eg. (1171) on the limiting trajectory is called the Lyapunov exponent, and a positive 
Lyapunov exponent indicates that the dynamical system is chaotic l|otl 120021) . The La grangian-history 



average of (Am) q " is different from the Lyapunov exponent in the choice of the director, but can be used 
to characterize potential ability of line stretching by the flow. From the viewpoint of distributive mixing, 
the magnitude of Am and its distribution in time and space can be used to characterize the mixing ability. 

While the line-stretching rate characterizes the stretching in one direction, the strain rate tensor 
describes the tri-axial straining states. By using the maximum line-stretching rate and D, the types of 
flow, namely the planar shear flow, uniaxial flow and biaxial elongational flows, are identified. The ratio 
of the maximum line-stretching rate to the s econd momen t of D is introduced in Eq. (|15[) . It is noted 
that fjj coincides with the mixing efficiency (Ottino, 19891) except for the choice of the direction n. By 
definition, fjj is an invariant of D. For incompressible flow, D : D is expressed by its largest and smallest 
eigenvalues, Am and A m , 



D 



D — \ M + X m 



+ (Am + A m y 



(19) 



therefore fo is found to be a function of A to /Am as 



\X M J 



(20) 



The incompressibility condition indicates that Am > 0, A m < 0, and the range of A to /Am € [—2, —1/2]. 
Thus fo is increasing monotonically on A to /Aa/ and lies in 



\<fn< 



1^ ^t<^ 



(21) 



For uniaxial and biaxial elongational flows, A,„/Am = —1/2 and A„j/Am = —2, respectively. For planar 
shear flow, X m /^M = — 1- The corresponding fjj is summarized in Eg. (|16p . 

Because the uniaxial elongational flo w establishes a more effective flow pattern in dispersive mixing 
than simple shear flow, several au t hors ([Yang and Manas-Zloczower , 1992 : Yao and Manas-Zloczower , 
Il998l : ICheng and Manas-Zloczower , 1997L Il998h have utilized some indicator of the strain pattern based 
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on D and Q like 



A = ^ , (22) 

2(£>:£>) 2 

= — 1 ( 23 ) 

o 

where f2 is the vorticity tensor and (.) indicates the steady co-rotational time derivative. Both A and Sf 
takes certain values for elongational flow, simple shear flow, and simple rotational flow. The flows can be 
discriminated based on the ratio of the a magnitude of vorticity $1 to a magnitude of V« which indirectly 
identifies the strain pattern: null vorticity, fi = 0, in pure elongational flow, and D : D = J7 T : fl in 
simple shear flow. By definition, A and Sf cannot distinguish between uniaxial and biaxial elongational 
flows. Like A and Sf, fn can be used to identify a type of strain-rate. It is noted that fo is different 
from A and Sf in that fo is based solely on the strain rate D and provides a direct characterization of D 
because we are interested in the stretching rate of the distance of two nearby points and the rotational 
component of flow does not affect the maximum line-stretching rate. 
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